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We consider a Fermi gas confined by a harmonic trapping potential and we highlight the role of 
the Fermi-Dirac statistics by studying frequency and damping of collective oscillations of quadrupole 
type in the framework of the quantum Boltzmann equation, in which statistical corrections are taken 
into account in the collisional integral. We are able to describe the crossover from the collisionless 
regime to the hydrodynamic one by introducing a temperature-dependent relaxation time tq. We 
show that, in the degenerate regime, the relaxation rate Tq 1 exhibits a temperature dependence 
different from the collision rate 7. We finally compare the collisional properties of the Fermi gas 
with the ones of the Bose gas for temperatures above the Bose-Einstein condensation. 
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C/3 ■ I. INTRODUCTION 

Since the experimental realization of Bose-Einstein condensation in trapped metastable atomic gases [0 the in- 
vestigation of collective oscillations has proven to be an important tool to characterize the behaviour of these many 
■ body systems. Indeed, the high accuracy of frequency measurements J2| and the sensitivity of collective phenomena 
to interaction effects (see || and references therein for a recent discussion of collective excitations in condensed Bose 
gases at finite temperature) makes them good candidates to check the predictions of theory Q . The quantum degen- 
' erate regime has also been recently reached in a magnetically trapped Fermi gas [^| and measurements of collective 
\Q , oscillations might be available soon also for fermionic species. 

The aim of this paper is to discuss the temperature dependence of the frequency and damping of the low energy 
collective oscillations of a trapped Fermi gas in its normal, non superfluid phase. The transition temperature to 
the superfluid phase is predicted to occur at very low temperature |^|,f^], while the effects of quantum degeneracy 
should be visible before the transition p|-|ri|]. For a discussion of the oscillations in the superfluid phase see |j7[[L2|. 
Collective oscillations can occur in the collisionless {lotq ^S> 1) or in the hydrodynamic (lutq <C 1) regime. Here cu 
is the frequency of the oscillation and tq is a characteristic relaxation time. We calculate the relaxation time tq by 
. linearizing the Boltzmann equation and we show that it is important to use tq instead of the inverse of the collision 
rate 7, which measures the time between collisions. The two times are proportional in the classical regime, but exhibit 
drastically different T dependences in conditions of quantum degeneracy. In addition to the study of the temperature 

O ■ dependence of tq in the Fermi case, which depicts the role of quantum statistics, we also compare the relaxation time 
• • 1 

of the collective oscillations with the one exhibited by a Bose gas above the critical temperature T c for Bose-Einstein 
condensation, to highlight the role of the bosonic statistics. 

The Pauli exclusion principle dramatically changes the collisional properties of Fermi gases in the quantum degen- 
erate regime and the onset of degeneracy can be traced down by looking at the corresponding changes of collective 
oscillations. At very low temperature Fermi-Dirac statistics quenches the number of collisions between atoms and one 
expects to find the gas in the collisionless (CL) regime. At higher temperatures there is the possibility of finding the 
gas in the hydrodynamic (HD) regime, especially if the cross section is very large, as in the case of 6 Li. An example 
of such a crossover has been given in |l]J for the spin dipole oscillation. 

The Fermi gas is taken to be a two component ("spin up" and "spin down") gas, to avoid the absence of s-wave 
collisions in the sample. The Bose gas will be instead considered to be trapped in a single spin state. We use a method 
recently developed by Guery-Odelin et al. in |l4| ] to evaluate in relaxation time approximation the collisional integral. 
The relaxation time tq is obtained through a linearization of the Boltzmann equation evaluated with an ansatz which 
describes the deformations in momentum space driven by the external oscillation. A similar technique has also been 
developed by Al Khawaja et al. |l5| for the evaluation of the collective excitations of a non-condensed Bose gas. An 
analytical and numerical investigation of the relaxation time and of the collision rate in a non-condensed uniform 
Bose gas was done by Lopez- Arias and Smerzi in |Tfj| ] while Kavoulakis et al. ]l7t presented a microscopic calculation 
of the relaxation time of a trapped non-condensed Bose gas based on the Boltzmann equation. Our approaches lead, 
in fact, to the same expression for the relaxation time. 

The paper is organized as follows: in Section || we briefly describe the method of the averages for a two component 
Fermi gas and we introduce the relaxation time tq for the quadrupolar collective oscillation and compare it with 
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the collisional rate 7. Section [I] is dedicated to the comparison with bosons. Section |y| reports the results for the 
collective oscillations of the Fermi gas as a function of the temperature. Section fv| draws conclusions and perspectives. 

II. RELAXATION TIME AND COLLISION RATE 

The equations describing the dynamical evolution of the quadrupole oscillations in a two component Fermi system 
can be derived starting from the Boltzmann equation 

'' Pl V + F a -V Pl )f a (r, Pl ,t)=I co u[f a } (1) 
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for the phase space distribution function for a given spin component a, f a {r, p, t) |18[] where F CT = — VV^ t (r) is the 
force originating from the external trapping potential, that we choose to be axially symmetric and that is, in general, 
different for different spin components 



^ t (r) = ^m[< ± (^+ y 2 )+<^ 2 ]. (2) 



The collisional integral reads 



l coll 



Ann^m J 

[(1 - /„(pi))(l - /-<r(p2))/ ff (pi)/- ff (Pa) - /«r(pi)/-«r(p 2 )(l - /a(pi))(l - /-<r(p£))] • (3) 

We have used the notation f a (Pi) = fa-(^,Pi,t). The moments p^ and p' 2 are constrained by energy and momentum 
conservation: p\ + p\ = p[ 2 + p' 2 2 and pi + P2 = p'i + P2- As the gas is at very low temperature the collisional cross 
section is taken to be c/ = Ana 2 with a the s-wave scattering length for up-down collisions. These collisions are 
the dominant ones at low temperatures, as up-up and down-down collisions occur in the p channel for antisymmetry 
requirements 19 1. By neglecting p-wave scattering, suppressed by the centrifugal barrier at the low energies of interest 
pOfl , we allow interactions only between up and down particles. 

The Boltzmann equation (Q) can be used to describe the trapped gas if one can use the semiclassical approximation 
and if mean field contributions to the energy of the gas are negligible. The former condition reads fc^T ^> Tiu)± tZ and 
can be fulfilled over a wide range of temperatures for large number of trapped particles as the Fermi temperature, 
below which quantum degeneracy effects begin to show, is defined as 7> = (SNy^hwho/kB (we have introduced the 
geometrical average of the trapping frequencies LJho — ( UJ '± UJ z) 1 ^)- The condition for the semiclassical approximation 
is then T/Tp ^> (SN)^ 1 ^ and is ensured if N is sufficiently large. Mean field effects are small the case for fermions 
p3| and this is also true for the Bose gas above the critical temperature for Bose-Einstein condensation T c = 
pV/£(3)] 1 ' 3 7ia;/, /fe.B. The reason is that the gas is dilute and the mean field interactions, of the order of gn(0) (here 
g = Ann a / '(2m,) with n(0) the central density), must be compared to the kinetic energy which, above T c , is of the 
order of ksT. The correction is typically of the order of few percent. 

When working with a two component system it is convenient to use the total distribution function / = /| + /j 
and the difference /_ = f] — fi- By summing and subtracting the Boltzmann equation for the up and the down 
component one gets 



[dt + — ■ V + F+ • V P1 ) /(r, pi, t) + F_ • V Pl /_(r, Pl , t) = I coU [ft] + hall [/J (4) 



in 



and 



'A + ?T ' V + F + ' V P! ) /- ( r > Pl > *) + F - ' V P! /( r > Pl . = T ooll [ft] ~ IcoU [ft] (5) 



in 



where F± = (F| ± FjJ/2. In the right hand side of equations (|J)-((|) the collisional contributions should be expressed 
in terms of / and /_ . In the simplified case of equal trapping potentials Vj xt = V} xt = V ext for the two spin components 
one has F_ = 0. If also iVj = N± then the difference /_ is zero, at equilibrium, for symmetry. In this work we will 
restrict our attention only to in-phase oscillations. One is then left with one equation, in which the only unknown is 
the total phase space distribution function / = 2/f = 2/|. 

A systematic way to solve the Boltzmann equation is provided by the Chapman-Enskog expansion [ pl| which 
allows to introduce the viscosity and the thermal conductivity in the hydrodynamic equations that will then exhibit 
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damping. To actually calculate the viscosity and the thermal conductivity it is necessary to provide some ansatz for 
their functional dependence. 

Starting from equation (^) with /_ = we instead derive equations for the average of generic dynamical quantities 
X(r,p): 

dt(x) ~ m-^p • V X ) - (F • V p x) = ]T(xW/<x]> (6) 

(7 

where the average is taken both in position and momentum space 

(x) = I X(r,p)/(r,p)dr I J /(r,p)dT (7) 
having defined dT = d 3 rd 3 p/h 3 . The collisional contribution can be rewritten in the form 

E<xWM) = iS / Ax(r,p)J co „ [/ CT (r, P )]dr / f f(r,p)dT (8) 

where Ax = xi + Xi ~ Xv ~ Xi' with x% = x( r :P«)- If X is a quantity conserved during the elastic collision then 
integral is zero. This happens if x = a ( r ) + M r ) • p + c(r)p 2 with a, b and c arbitrary functions of the position. 

As shown in [p"i| the dynamical evolution of the m z = quadrupole oscillation in an axially symmetric trapping 
potential is described by a system of six coupled equations that we rewrite here: 

d t (xi) ~ 2<X3> = 
d t (x 2 )-2(x 4 ) =0 

ft<X3> - <X5> + ^f^(xi) + ^^ife) = 

te> - (Xe) + 3 ^ (X0 + ^^(X2) = (9) 

_ . . 2iv 2 z + 4u) 2 ,. . 2u 2 z - n 
dt(X5) + ~(X3) + -(Xi) = 

a t <Xe) + ^(xs) + ^<X 4 ) = 2^<Xe/coii[/«r]). 

The Xi quantities are defined as 

2 r, 2 2 

Xi = r X2 = 2z - r ± 

X3 = r • v X4 = 2zu 2 - r_L • v_l (10) 

2 r, 2 2 

X5 = U X6 = 2w z - U ± 

where v = p/m and rx, vi are respectively the projections of the r and v vectors in the xy plane. The m 2 = 
quadrupole oscillation is coupled to the monopole one unless the trap is spherical (oj z = u>±). The X6 combination is 
the only one for which the collisional contribution does not vanishes. To describe the m z = 2 quadrupole oscillation 
one should instead introduce the mean values oi xy, v x y + v y x and v x v y . We shall not consider here this mode. 

Solving system (||) is still a hard task because a full solution of the Boltzmann equation is needed to actually 
calculate the average of the collisional contribution. What can be done is to linearize the distribution function 
assuming a specific form for its deviation from equilibrium, finding a variational estimate for /(r, p). By recalling 
that in the collisionless regime one can have solutions with anisotropic equilibrium velocities, corresponding to different 
temperatures in the radial and axial directions, we make the ansatz 

/ CT (r, p) = / (r, (1 + o^) 1 / V, (1 + «,) 1/2 P,) (11) 

where /o is the equilibrium distribution function. The proposed anisotropy in momentum space is of quadrupole type 
if one chooses «i = —a z /2. In this case the normalization is preserved to first order in the deformation parameters. 
Deformations in the real space do not need to be taken into account in ansatz (|ll| ) because, after linearization, they do 
not contribute to the integral as a consequence of the invariance property of the collisional integral already discussed. 
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The linearization of f a around a±_ = a z = leads to 

/^/o-(^ + ^V-/o) (12) 
V 2m 2m J 

with (3 = (ksT)- 1 and, using the equilibrium identity (1 - /o(pi))(l - /o(P2))/o(p'i)/o(p 2 ) = /o(pi)/o(P2)(l - 
/o(p'i))(l — /o(P2))j we nn d the following expression, to first order in the alpha's, for the collisional term: 



E 



(Xelcoidfa]) - 16 %f mN J d 3 rd 3 Pl d 3 p 2 <mA Xf >\vi-v 2 \ x 



(a ± Ap 2 ± + a z Ap 2 z ) / (pi)/o(p 2 )(l - /o(p'i))(l - /o(p 2 ))- (13) 

Carrying out the calculation (see Appendix) shows that the collisional contribution is proportional to the mean value 
of (xg) so that one can write 

1 Z{X»Icoll\f<r]) = -— (14) 

defining the temperature dependent relaxation time for the quadrupole oscillation. The dimensionless quantity UhoTQ 
can be written in the form 

4 V=^)Vr) us. 



UhoTQ 5\/3n V a hoJ \Tp J 

where a/ lo = [%/ (mmho)] 1 ^ 2 is the harmonic oscillator length and Fq{£) (see Appendix) is an universal function of 
the reduced temperature t = T/Tp, showed in Figure |l|, which contains the effect of Fermi statistics. Numerical 
coefficients have been chosen so that Fq has the limiting form 1/t for t ^> 1 (classical regime). In the fully degenerate 
regime the integral defining Fq is analytical and the result 

F Q (t) = 8iT 2 t 2 (16) 

is found. The T~ 2 behaviour for the relaxation time is typical of degenerate Fermi systems |2^j23|| . We recall that 
result ( |l6| ) holds above the BCS phase transition. Conversely we point out explicitly that result ( |15| ) in the classical 
limit is correct only as long as the approximations that led to the collisional integral (||) remain valid. The truly 
classical system is the one for which all partial waves are taken into account in the two body collision, thereby 
restoring the collisions between up-up and down-down particles, no longer suppressed by (quantum) antisymmetry 
requirements. The classical limit of equation ( |l5| ) instead keeps on having the quantal characteristic of up-down 
s-wave scattering only. However p-wave collisions in real experiments become important only at temperatures much 
higher than the Fermi temperature |pfj|| and it is safe to describe in such a way our system. We finally observe that 
the relaxation time depends on the square of the scattering length a and, consequently, it is not sensitive to its sign. 
Mean field corrections, that have been neglected in the Boltzmann equation, do depend on the sign of the interaction 

II- 

Using result ( |l4| ) the system of equations (g) is closed and the dispersion law for the collective oscillations can be 
obtained. Before considering the oscillation frequencies we address our attention to the number of collisions that take 
place in the gas by calculating the collision rate 
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J d 3 rd 3 Pl d 3 p 2 dn |vi - v 2 |/ 1T / 2i (l - A, T )(1 - f vl ). (17) 

The factor 2 in equation ([l7]) takes into account the presence of the two components. By treating integral ( |l7| ) in the 
very same way of the collisional contribution ([l5]) we find the following expression for the rate at equilibrium (put in 
adimensional form): 



-2- = — !— = — L- NV 1/3 — ) F 7 f — 

Who ^hoTj v37T V a ho J \Tf 



(18) 



where again the coefficients of the function F 1 are chosen so that it has the 1/t classical behaviour. In the degenerate 
regime one instead finds the analytical result 

' 3 2 / 3 8tt(V/3 J-) 2 (jL)\ (19) 



Who?-) t->o \ aho J \Tf J 
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Result (|l9|) explicitly shows the importance of using the relaxation time tq instead of the inverse of the collision rate 
t 7 : while in the classical regime they are proportional, their ratio being (rQ/r 7 ) c ; ass . = 5/4, in the degenerate regime 
their temperature dependence is different. The cubic dependence on temperature of the collision rate in the quantum 
degenerate regime is not a consequence of the trapping but holds also in the uniform system where the same type of 
analytical calculations can be carried on as well. 



III. COMPARISON WITH THE BOSE CASE 



The absence of proportionality between collision rate and relaxation time in a non condensed Bose gas was pointed 
out by Lopez-Arias and Smerzi |l6| who numerically calculated the relaxation time tq for the uniform gas. The 
method of averages allows to calculate tq for the trapped Bose gas by considering the same ansatz (|ll|) used in the 
Fermi case. We now have only one spin component and thus the spin index is absent. The cross section is erf, = 87ra 2 
because the scattering is between indistinguishable particles and the collisional integral reads 

/co " [/] = 4^fc J d 3 p 2 dn\ Pl -p 2 \ x 

[(1 + /( Pl ))(l + /(p 2 ))/(pi)/(pi) - /(Pi)/( P2 )(1 + /(Pi))(l + f(p' 2 ))\ ■ (20) 

The (1 + /) factors in the above equation account for the increase of probability of a scattering event in an already 
occupied final state, typical of Bose systems. 

The expression for the quadrupolar relaxation time Tq , analog to ( |l5|) and in agreement with |T?j| , is given by 



16 VC(3) 1/3 o V^i ( T 



UhoTn 5-7T I a ho J v \T C 



Q 

and the one for the collision rate, analog to (|lj), reads instead 



^ 1/3 — ) G Q Nr ( 21 ) 



76 _ 1 _4j/g3) f Ml/3 a \* (T 



N^— G, b l±- . (22) 



Uho u ho r lb IT \ ahoJ \J-C 

The functions Gq and G lb are plotted in Figure || and their expressions are given in Appendix. Both functions behave 
as T c /T in the classical limit. Comparison of the two equations shows again the limiting ratio (tq/t 76 ) c / qss . = 5/4 
fl4| . By using an improved trial function the authors of jnj have shown that in the classical limit the ratio changes 
by only 1%, proving the validity of ansatz (Ej 



The quantum behaviour of Tq and r 76 is opposite to the behaviour of the corresponding fermionic quantities. 
Quantum statistical effects enhance the number of collisions taking place in the trapped gas and the rate 7& is seen 
to increase with respect to the classical trend. The integral G lb of equation ( p2| ) (see Appendix) is well behaved as 
Pfj, — > 0~ and reaches a finite value, differently from the uniform case Jl6[], where the collision rate diverges at the 
critical temperature, driving to zero the time between collisions. This divergence has also been found by Nikuni et 
al. ]2~i[ ], in the framework of a quantum Boltzmann equation coupled with an equation for the condensate, as T c is 
approached both from above and from below. In the non-uniform case instead r 7t does not vanishes at the transition. 
To understand this it is useful to recall that r 76 , as well as Tq and their fermionic counterparts, can be regarded as 
the average over the whole sample of a local collision (relaxation) time. It is then natural to find a weakening of the 
divergences exhibited by uniform systems. Figure || shows that the quadrupolar relaxation rate (tq) -1 , similarly to 
76, increases its value with respect to the classical value but in a less stronger way as compared to the increase of the 
collision rate. This witnesses the loss of efficiency of the collisions near the transition temperature: the number of 
collisions gets larger, due to the Bose factors (1 + /), but the relaxation rate does not increases its value proportionally 
p5| . This feature was pointed out also in [flql where the authors even find a diverging behaviour for the relaxation 
time of the quadrupolar oscillation. To obtain a diverging behaviour in our trapped case one should observe the 
vanishing of the function Gq in equation (|l]) as T — ► T c + . This is not the case, as the integral is positive defined 
and increases its value near the critical temperature. As final remark we observe that the (linearized) Boltzmann 
equation is not accurate in describing the system at the Bose-Einstein phase transition because it completely neglects 
fluctuations and coherence phenomena so that one should not look at the quantum limit (3/i = as a quantitative 
prediction. Conversely the behaviour of Gq and G lb as T — > T c + is physically relevant. The coupled dynamics of the 
condensate and the thermal component in the trapped case below T c can be found in [p6| where the normal component 
and the Bose condensed fraction are taken into account in the framework of a quantum Boltzmann equation coupled 
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with a generalized Gross-Pitaevskii equation. Finally two examples of theories that also deal with the issue of the 
condensate formation, within formalism derived from quantum optics and from field theory, can be found in |27|] and 




IV. QUADRUPOLE OSCILLATIONS IN THE FERMI GAS 

We now consider the frequency and the damping of the quadrupolar collective oscillation that can be derived from 
the set of equations (^|) together with ([u]) . By looking for solution of oscillatory type one finds the dispersion relation 
by imposing the determinant of the system to be zero. The dispersion relation coincides with the one already derived 
for the classical gas Jl4[ , as quantum statistical effects are kept into account in the relaxation time tq . One gets 

.2 A, ,2\ I, .2 a. ,2\ 1 (, A 2 ,.2/V. ,2 i a. ,2\ , o, ,2 , .2 



(uj z - Auo z z ) - 4cj1) - — I cj 4 - -uj z (5cj1 + Auoi) + 8^1^ J = 0. (23) 

The first term of (E^) corresponds to the dispersion law for the pure collisionless regime (lutq — > 00). In this 
case the eigenfrequencies coincide with the ones predicted by the single particle harmonic oscillator Hamiltonian: 
locl = 2wi.^. The total Hamiltonian is in fact the sum of single particle Hamiltonians and interactions are completely 
neglected. In this case the concept of collective oscillation actually looses his meaning. On the other hand the term 
multiplying {/(lotq) refers to the pure hydrodynamical regime (lotq — > 0) and the dispersion relation gives the result 
p9| uJjjjj/uj 2 ^ = (4A 2 + 5 ± \/l6A 4 — 32A 2 + 25)/3 where A = uj z /loj_ is the anisotropy parameter. Putting A = 1 gives 
two frequencies uj^ d = 2u>ho and oj^ d = \f 7 iui}- lo corresponding to the monopole and quadrupole mode respectively 
(as observed, system (|9|) decouples in two sub-sets). The relevant limits for deformed traps A <C 1 and A 3> 1 
can be obtained as well Jl|,[l^] • When the relaxation time tq is of the order of the oscillation frequency one is in 
an intermediate regime and equation (p3|) is solved numerically for u> — to r + iT yielding the frequency uj t and the 



damping T of the oscillation. As observed in 30 in the Bose gas case, the notion of hydrodynamic regime in trapped 
gases fails at the boundary of the trap, where the gas is rarefied and collisions are rare. The average performed over 
all the sample in calculating the relaxation time tq takes care of this effect. Indeed the possibility of an unified 
description of the collisionless and of the hydrodynamic regimes is the useful feature of the presented model. 

By looking at the temperature dependence of the relaxation rate, depicted in Figure [jj it is clear that the gas will 
be in the collisionless regime both at low and high temperatures. At low temperatures the Pauli principle reduces the 
number of available final states for the scattering process and the relaxation rate tends to zero as (T/Tp) 2 . At high 
temperatures, as noted previously, the relaxation rate is proportional to the collision rate and decreases as Tp/T. 
This is basically a density effect: for fixed number of particles a Boltzmann gas gets more and more rarefied as 
temperature increases. In the intermediate region there is room for an increase of Tq 1 . Whether the relaxation time 
tq will reach a sufficiently small value to bring the system in the hydrodynamic regime depends on the combination 
(A 1 / 3 a/a/ lo ) 2 of the number of particles, the scattering length and the trap parameters of equation (|l5| ) as the function 
Fq depends on the reduced temperature only. The maximum of Fq is ~ 1.8 at T/Tp ~ 0.35 so that a minimum 
value of uJhoTQ = 3.15/(N 1 / 3 a/aho) 2 is reached. Large numbers of trapped atoms as well as large scattering lengths 
should allow the observation of the hydrodynamic regime. In the following we will consider explicitly the relevant 
experimental cases of 40 K || and 6 Li 

The crossover from collisionless to the hydrodynamic regime is most striking in the case of the spin dipole oscillation 
fl3"f where the dispersion law is simpler than ( |23| ) and allows to discriminate analytically between a damped oscillatory 
behaviour and an overdamped one. In the case of spin dipole oscillations one can find overdamping because the spin 
current Vf — v_[ is not conserved during collisions and this leads to a purely diffusive mode in the collision dominated 
regime. In the case of quadrupole oscillations instead both the collisionless and the hydrodynamic regimes support 
an oscillation and the actual regime can be deduced from the value of the real part of the oscillation frequency which 
passes, e.g. in the spherical case, from 2ojho to y/2u>ha- The spectrum of a trapped Fermi gas in the hydrodynamic 
regime has been studied independently by Amoruso et al. |}2| and by Bruun and Clark |3^] for T <Z^Tp and T ^> Tp. 
The frequency of the quadrupolar oscillation clearly coincides with result (^3j) in the lotq — ► limit. However 
we observe that in the limits of low and high temperature the value of the oscillation frequency will not be the 
hydrodynamic but the collisionless one. This even supplies a possible signature for the superfluid transition since the 
value of the quadrupole frequency for a superfluid Fermigas coincides with the hydrodynamic value [Q and at the 
calculated temperature for the BCS phase transition ||(7| tq is large and the system is in the collisionless regime. 
So for large number of particles and for large values of \a\/a,ho one should observe, as a function of the decreasing 
temperature, the sequence of collisionless - hydrodynamic - collisionless - hydrodynamic value for the real part of 
the quadrupole frequency and a corresponding structure of local maxima and minima in the damping rate of the 
oscillation. 
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To summarize this discussion we show in Figure |3j the real and the imaginary part of the quadrupolar oscillation 
of a spherical (A = 1) cloud of 40 K atoms. The scattering length is [0 a = 157 ao with do the Bohr radius. We take 
the ratio a/a,ho ~ 10 -3 and N = 10 6 . The minimum value of LOhoTQ is then ~ 350 and the system is always in the 
collisionless regime. In the upper panel the real part of the frequency (solid line) is compared with the frequency of 
a Boltzmann gas (dashed line) which, at low temperatures, is always in the hydrodynamic regime. This emphasizes 
the role of Fermi statistics in the low temperature regime of the trapped gas. The lower panel shows the damping 
of the oscillation. The two gases exhibit zero damping in the pure collisionless regime at high temperature. At zero 
temperature both the Boltzmann gas and the Fermi gas have zero damping but the former is in the pure hydrodynamic 
regime while the latter is in the pure collisionless one because of Pauli blocking. Inspection of the damping rate gives 
then another clear signature of quantum degeneracy as it departs significantly from the classical value. 

Figure |] shows instead the real (upper panel) and the imaginary (lower panel) part of the quadrupolar oscillation 
frequency for N = 10 6 atoms of 6 Li, that have a resonant value of the scattering length a = —2160 ao, in a spherical 
trap. Due to the lighter mass with respect to potassium the ratio \a\/ah is ~ 5 x 10~ 2 and a minimum value of 
WhoTQ ~ 0.14 is found. Lithium can thus be found in the hydrodynamic regime in the region around the maximum 
of the function Fq (Figure [l]). Indeed the real part of the frequency bends towards the hydrodynamic value \f2u>h 
as long as the Fermi statistics allows a sufficient collisional activity. Correspondingly the damping has a relative 
minimum. One should point out that if the system had reached the pure hydrodynamical regime at T ~ 0.35 Tp the 
damping would have then been exactly zero. As the temperature is further lowered Pauli principle forces the system 
to become again collisionless, the real part raises again towards the value 2 u>ho while the damping exhibit a local 
maximum and then falls to zero once the complete collisionless regime is reached. 

The two proposed examples involve the measured scattering lengths for potassium and lithium. Different predictions 
for the frequency and the damping can be obtained for different values of the number of trapped atoms, as well as of the 
trapping parameters and scattering length through the use of the definition (|l5| ) of the relaxation time together with 
the universal function plotted in Figure |l|. A way to reach the hydrodynamic regime also in the case of the potassium 
isotope could be the change of the value of its scattering length with external fields, as already experimentally observed 
in alkali atoms pi ]. 



We have studied the collective frequencies and the damping of a trapped normal Fermi system through the use of 
a linearized Boltzmann equation, a method that gives a variational estimate for the relaxation time tq relative to 
the coupled monopole-quadrupole oscillations. The effects of Fermi statistics are included in the relaxation time tq 
that diverges at zero temperature as a consequence of the Pauli principle and at classical temperatures T > Tp as a 
consequence of the diluitness of the gas. In both cases the collisionless regime for the collective oscillation is reached. 
At intermediate temperatures we have found the possibility of a hydrodynamical behaviour of the quadrupole mode for 
the 6 Li isotope, due to his large scattering length, for reasonable number of trapped particles. We have also compared 
the relaxation rate Tq 1 with the collision rate 7, finding analytically different behaviours at low temperatures. This 
emphasizes the need of distinguishing tq from 7 _1 . The role of quantum statistics has also been considered in the 
case of a Bose gas and the difference between the collision rate and the relaxation rate have been pointed out also 
in this case. Future developments will concentrate on non symmetric configurations, that may be relevant from the 
experimental point of view, as well as on the study of the damping of collective oscillations in mixed systems. 

Useful discussions with A. Smerzi and H. T. C. Stoof are acknowledged. I wish to thank S. Stringari for critically 
reading the manuscript. 



V. CONCLUSIONS 



APPENDIX A: CALCULATION OF THE COLLISIONAL INTEGRAL 



We first introduce the center of mass variables 



v 1= V + - 
v 

v 2 = V - ~ 

2 




(Al) 



v' 

v' = V - — 

2 2 



7 



where |v| = |v'| = v for momentum conservation. The Jacobian associated to transformation (Al) is one. In these 
variables one has 

1 



2 v?) 



1 



and Ax 6 = 2Av 2 — Auj_. The product /o(vi)/o(v 2 )(l — /o( v 'i))(l — /o( v 2)) can be written in the useful form 
/o(vi)/o(v 2 )(l-/ (vi))(l-/oK)) 



(A2) 



(A3) 



4 cosh(^) + cosh(7V • v) cosh(^) + cosh(7V • v') 
where v = (3(m(V 2 + i> 2 /4)/2 + mui 2 ± (x 2 + y 2 + X 2 z 2 )/2 — //) and 7 = f3m/2. We then have to linear order in a±, a z 

a 

with 5a = a z — a 1 and 



/ = y d 3 rd 3 Vd 3 vdnv (v 2 z -v' 2 Y 



cosh(zy) + cosh(7V • v) cosh(^) + cosh(7V • v') 



(A5) 



Integral (A5) can be carried on as follows: introduce a reference frame (a, b,c) oriented such that V = Vk and z 
is in the plane generated by a and b. We then have the following expressions: 



v z — v a{k ■ z) + Vb(b ■ z) = v cos($)(a • z) + t>sin($) cos((^)(b ■ z) 
v' z = v' a (k ■ z) + v' b (h ■ z) = v cos(i?')(a ■ z) + v sin(t/) cos(y/)(b • z) 



(A6) 
(A7) 



and V • vW = Vv cos(z9 < -'- 1 ). Integrations in dip and dip' are immediate. One can then integrate in dd and d"D' using 
the fact that cosh(a;) is an even function to end with the following intermediate form of the integral 



I = 



2V 
~5A 



■ J r 2 drV 2 dVv 7 dv J dxdy Q + | 



2 22 
x — x y 



1 



1 



cosh(y) + cosh^Vux) cosh(^) + cosh(7l/wy) 



(A8) 



A final change of variables to polar dimcnsionless coordinates leads to the final form of I which has to be integrated 
numerically except in the low and high temperature limits: 



I = 



2 21 7T 4 



1 



5co 3 ho (p m y Jo 



p 6 dp / dddip cos 2 sin 10 sin 7 (<p) cos 2 (ip) x 



/: 



1 



cosh(/9 — (3p) + cosh(psin (t?) sin(2^)?;) 

1 



(A9) 



The value of (x6) can be calculated to be, to first order in 5a 

, . 4(5a 
(Xe) 



Nm0*(hu ho ) 3 



cosh(/9 — + cosh(p sin 2 (•&) sin(2</j)£) 
U(z) (A10) 



and the common proportionality of (A5) and ( A10 ) on 5a allows to write equation (|14j) of the text. Using ( A10 ) and 
(A9) one builds the Fq function as 



F c 



T \ 



9 2 5 



TfJ h{z) \Tf 
1 



T 



p 6 dp I d^cos 2 (i?)sin 10 (i?)sin 7 (^)cos 2 (^) x 

1 



cosh(/9 — I3p) + cosh(psin (1?) sin(2ip)ri) 

1 



cosh(p — (3fi) + cosh(psin ($) sin(2y>)£) 



■ (AH) 



8 



In the classical limit [3p — > — oo and the denominators such as cosh((p — j3p) + cosh(psin 2 ($) sin(2</?)£) can be 
approximated with exp(p — [3p)/2. In the quantum degenerate regime instead ftp — > oo and, after changing variables 
to p' — (p — /3fi), f — (p' + f3p) sin 2 ($) sin(2<^)^ and 77' = (p' + (3p) sin 2 ($) sm(2(p)rj one has integrals over the whole 
real axis that can be done analytically. 

The same type of calculations can be done for the collision rate as well as for the corresponding bosonic expressions 
of Tq and r 76 . The formula for F 1 is: 

t \ a 9? / T 1 \ 5 r°° r, 2 W 2 W 2 



JO 



atanh 2 (tanh ((p - /3^)/2) tanh (psin 2 (i?) sin(2<p)/2)) . (A12) 



In the Bose case one finds the Bose function 34(2) in place of the Fermi function fi{z) and the expressions for the 
bosonic quadrupolar relaxation function Gq and for the bosonic rate function G 7b are 



° Q (j) = c(3W(*) (s) r p6dp i 2 ^ cos2 w sinio w sin7 ^) c ° 82 ^) x 



^2 2 /"2 



cosh(p — /3/i) — cosh(,o sin 2 (1?) sin(2<p)rf) 

o , (A13) 

cosh(p — ftp) — cosh(psin (■&) sin(2<£)£) 



G f^-) = -|— f—) 5 r ^ r /2 d$dtp sin 2 (20)cosfo>) x 

7b Vr c y ttC(3) 2 vTJ i sinh 2 (p-/3 M )7o ^ 

atanh 2 f tanh(psin 2 W sin(2 y )/2)y 

V tanh((p-/3/x)/2) J y ' 

The classical limit is performed again by taking (3p — > —00. Differently from the Fermi case, as the critical temperature 



T c is approached the fugacity tends to one, i.e. [3p — > 0. In this limit the integrand in (A14) diverges in the origin as 



a logarithm but the integral is finite. This is a consequence of th e confinement, since in the uniform case the collision 



rate diverges at the critical temperature [^6|j24j]. Integral ( A13 ) can be worked out analytically to include the effect 



of quantum degeneracy to first order in the fugacity z — exp((3p) far from th e tra nsition temperature (t 3> 1). We 
perform the expansion to check explicitly the equivalence of our expression ( A1S ) with the one presented in [fl7| . 



By expanding the denominator in the integral in powers of the fugacity and re-introducing Cartesian coordinates we 
obtain 



o q(t) 128tV ' 



Q{2>)irgi{z 
which gives 



/>oo />1 

/ x\dxzx\dxix\dx\ I d-qd^ (l + 2r] 2 — 3?7 2 £ 2 ) x 
{exp[-2(x 2 + x\ + xl)] + 2zexp[-3(.T 2 + x\ + a; 2 )] (cosh(2x 1 x 2 r]) + cosh^xi^O)} (A15) 



t 2 7 2 { 3 

G Q {t) ~ -1 . . 1 + ~* I . (A16) 



C(3) 34 (z) V 8 

From the normalization condition g^(z) = £(3)/£ 3 one obtains the correction to the fugacity 



which, together with the expansion of 34(2) — z + z 2 /16 lead to 



1 A 3 C(3) 



G °®*it 1+ 16^) (A18) 







which is the classical limit 1/t with the quantum correction of |M . A similar expansion can be performed for the rate 
function G lb as well as for the functions Fq and F 1 . We report the result for the Fq function which is 
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FIG. 2. Quadrupole universal function (solid line) and rate function (dashed line) for the Bose case as a function of the 
reduced temperature T/T c . The classical limit 1/t is shown as a dotted line. 
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FIG. 3. Real and imaginary part of the quadrupole oscillation frequency for N = 10 6 40 K atoms (solid lines). The dashed 
lines are the corresponding values for a classical gas. 
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